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Abstract 

The first order nonlinear ODE ip(t) + sin ip(t) = B + A cos ut, (A, B, oj 
are real constants) which is commonly used as a simple model of an 
overdamped Josephson junction in superconductors is investigated. 
Its general solution is obtained in the case of the choice of parameters 
associated with one of three possible kinds of asymptotic behavior of 
solutions known as phase-lock where all but one solutions converge 
to a common 'essentially periodic' attractor. The general solution 
is represented in explicit form in terms of the Floquet solution of a 
particular instance of the double confluent Heun equation (DCHE). 
In turn, the solution of DCHE is represented through the Laurent 
series which defines an analytic function on the Riemann sphere with 
punctured poles. The Laurent series coefficients are given in explicit 
form in terms of infinite products of 2 x 2 matrices with a single 
zero element. The closed form of the phase- lock condition is obtained 
which is represented as the condition of existence of a real root of 
the transcendental function. The phase-lock criterion is conjectured 
whose plausibility is confirmed in numerical tests. 

1 Introduction 

The nonlinear first order ODE 



<p(t)+ sin tp(t)=q(t), (1) 

is commonly used in applied physics as the simple mathematical model de- 
scribing the electric properties of Josephson junctions (JJ's) in superconduc- 
tors [H Here the RHS function q(t) assumed to be known specifies the 
external impact to JJ representing the appropriately normalized bias current 
(or simply bias) supplied by an external current source. The unknown real 
valued function ip(t) called the phase describes the macroscopic quantum 
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state of JJ. In particular, it is connected with the instantaneous voltage V 
applied across JJ in accord with the equation V = (h/2e)d(p/dr, where h~ is 
the Plank constant, e is the electron charge, r is the (dimensional) current 
time. The dimensionless variable t entering Eq. ([1]) is defined as t = uj c t, 
where u c is a constant parameter depending on the junction properties and 
named JJ characteristic frequency. See Ref [5] for more details of JJ physics. 

Eq. ([I]) arises as the limiting case of the second order ODE utilized in more 
general Resistively Shunted Junction (RSJ) model [31 Hj. The reduction is 
legitimate if the role of the junction capacitance proves negligible. In practice, 
if JJ can be described by Eq. ([1]) it is named overdamped. Summarizing the 
aforementioned relationships, we shall name Eq. ([T]), for brevity, overdamped 
Josephson junction equation (OJJE). 

Under concordant conditions, the theoretical modeling applying OJJE is 
in excellent agreement with experiments. It is also worth noting that nowa- 
days electronic devices based on the Josephson effect play the important 
role in various branches of measurement technology. In particular, JJ ar- 
rays serve the heart element of the modern DC voltage standards [5J. The 
development of JJ-based synthesizers of AC voltage waveforms is currently 
in progress [3 [HJ [9] . These and other successful applications stimulate the 
growing interest to the theoretical study and the modeling of JJ properties 
including investigation of capabilities of RSJ model and its limiting cases and 
the predictions they lead to. 

To be more specific, the case most important from viewpoint of applica- 
tions and simultaneously distinguished by the wealthiness of the underlaid 
mathematics is definitely the one of the bias function representing harmonic 
oscillations. Without loss of generality, it can be get in the following form 

q(t) = B + Acosujt, (2) 

where A,B,u > are some real constant parameters. Hereinafter, the ab- 
breviation OJJE introduced above will refer to the couple ([I]),©. 

In spite of apparent simplicity of OJJE, few facts of its specific analytic 
theory had been available until recently. Some preliminary results concerning 
the problem of derivation of analytic solutions of OJJE in general setting had 
been obtained in Ref. [12]. The approach put forward therein is elaborated in 
the present work. The focus is made on the case of manifestation of the phase- 
lock property by Eq. ([T]) which is one of its most important features from 
viewpoint of applications. The phase-lock property is formalized as follows: 
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in the case of phase-lock any solution tp = ip(t) of Eq. ([T]) either yields 
a periodic exponent e 19 , exactly two such distinguished solutions existing, 
or, as the time parameter grows, e llf exponentially converges to the similar 
exponent for the one (common for all <£>'s) of the periodic functions just 
noted (another one is the repeller). (Here and in what follows we shall not 
distinguish phase functions which differ by a constant equal to 2n times an 
integer.) The corresponding period coincides with one of q(t), i.e., in the 
case (T5]), 27rc<j~ 1 . This behavior is stable with respect to weak parameter 
perturbations, i.e. the subset of parameter values leading to phase-lock is 
open. It is worth noting for completeness that in the opposite (no phase- 
lock) case no stable periodicity in the behavior of e ltp is observed. There is 
also a third, intermediate type of the phase behavior, where the attractor 
and repeller are, in a sense, merged. It is realized on the lower-dimensional 
subset of the space of parameter values [11] . 

In the present work, the complete analytic solution of OJJE is obtained 
under assumption of the parameter choice ensuring the phase-lock property. 
The closed form of the phase-lock criterion in the form of the constraint 
imposed on the problem parameters is conjectured. 



2 Overdamped Josephson junction equation 
against reduced double confluent Heun equa- 
tion 

The analytic theory of OJJE ([I]),© can be based on its reduction to the 
following system of two linear first order ODEs [T2] 

Aiuz 2 x'{z) = 2zx{z) + [2Bz + A(z 2 + l)]y{z), 
-Awz 2 y\z) = [2Bz + A(z 2 + l)]x(z) +2zy(z), (3) 

where z is the free complex variable. Indeed, on the universal covering Q\ ~ 
M 3 t of the unit circle in C, i.e. for z = exp(iu;t), any non-trivial solution of 
Eqs. (ED determines a solution of OJJE in accordance with the equation 

— i^fty 

6XPM = WxTWy- (4) 

(3? denotes the real part) supplemented with the continuity requirement. 
[The fulfillment of OJJE follows from a straightforward computation taking 
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into account the equality 1 = z~ l holding true on Qi which is utilized for the 
demonstration that, for real A, B, u, the functions dix, dty also verify Eqs. (J3]) 
on Q\.] Conversely, any real valued solution (p(t) of OJJE induces through 
Eq. (j3J) some 'initial data' x(0), y(0) (with arbitrary norm (x 2 (0)+y 2 (0)) 1 ^ 2 > 
0) for Eqs. ([3]). Having solved the latter on Qi for x, y, one obtains, applying 
OH), another phase function obeying OJJE. It however must coincide, due to 
the identical initial value assumed at t — 0, with the original tp(t). Finally, 
since these x, y defined on Qi simultaneously obey the linear ODEs ([3]) with 
meromorphic coefficients which have the only singular points z = and 
z^ 1 = 0, they can be extended, integrating ([3]) along radial directions, to 
analytic functions defined on the whole Q = Qi x M + , the universal covering 
of the Riemann sphere with the punctured poles z = and z~ x = 0. 

It is worth noting that on Qi the real valued functions x = x(t) = 
SR x (e 1U3t ),y = y(t) = $t,y(e luJt ) satisfy the equations 

2 d x/dt = x + qy, —2dy/dt = qx + y. (5) 

leading, together with (J3J), to the equation (d/dt)log(x 2 + y 2 ) = cos<£> which 
implies 

consti e~* < \x + iy\ 2 < const2 e* for t > 0. (6) 

Notice that if x, y ^ then consti, const2 may be assumed to be strictly 
positive. We shall refer to these boundings later on. 

The key observation enabling one to radically simplify the problem of 
description of the space of solutions of Eqs. ([3]) is as follows [12]. Let us 
introduce the analytic function v(z) which satisfies the equation 



d^ , , „ s d 



^ + (M^ 2 + i)-^)^ + (2^)- 2 



0, (7) 



where 



are the constants replacing original A, B which will be used below whenever 
it proves convenient. Then a straightforward calculation shows that the 
functions x, y determined by the equations 

v = i z^~ exp f-fi (— z + -z^ 1 )^ (x — iy), (9) 

v' = (2u;z)~ 1 z Ii f L exp Q/i (-z + z^ 1 )^ (x + iy) (10) 
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verify Eqs. fl3]). Conversely, defining the function v(z) through the solution 
x, y of Eqs. ([3]) in accordance with Eq. flUD, a straightforward computation 
proves satisfaction of Eq. ffTOl) and, then, Eq. ([7j) follows. Thus, Eqs. (j3J) are 
equivalent to (jTJ) and Eqs. fl9l). f|T0l) represent the corresponding one-to-one 
transformation. 

Eq. (j7j) coincides, after appropriate identification of the constant param- 
eters, with Eq. (1.4.40) from Ref. [T5]. It represents therefore a particular 
instance of the double confluent Heun equation (DCHE) which can be shown 
to be in our case non-degenerated. 

It is also worth reproducing here the canonical form of the "generic" 
DCHE as it is given in Ref. [T5] (Eq. (4.5.1)). It reads 

z 2 pi + (-z 2 + cz + t)p- + {-az + X)y = (11) 
dz 2 dz 

where a, c, t, A are some constants. To adjust it to our case, a single term 
has to be eliminated setting a = 0. For brevity, we shall name this subclass 
of DCHE's reduced. Besides, some obvious rescaling of the free variable z 
is to be carried out. After these, the three constant parameters remained 
in the resulting equation correspond to our constant parameters n, fi, u (we 
shall not need and so omit the reproducing of the concrete form of this 
transformation) . 

The general analytic theory of DCHE is given in the chapter 8 of the trea- 
tise [H] . DCHE solutions are there represented, up to nonzero factors given 
in explicit form, through the Laurent series whose coefficients are assumed 
to be computable through the 'endless' chain of 3-term linear homogeneous 
equations ('recurrence relations'). In the present work, we derive the solution 
of reduced DCHE in a cognate but more explicit form. 

As a technical limitation, we also stipulate in the present work for the 
additional condition to be imposed on the free constant problem parameters 
claiming of them the ensuring of the phase-lock property. On the base of 
practice of numerical computations, it can be conjectured that such param- 
eter values fill up a non-empty open subset {phase-lock area) in the whole 
parameter space (see also the Conjecture A below). The case where the pa- 
rameters belong to its complement is left beyond the scope of the present 
work. 
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3 Formal solution of reduced DCHE by Lau- 
rent series 



Let us introduce yet another unknown function E(z) replacing v(z) by means 
of the transformation 

v ( z ) = z^- iie e-> iZ E(z), (12) 

where the discrete parity parameter e may assume one of the two values, 
either e = or e = 1, and x is some real positive constant which will be 
determined latter on. For v obeying (171) , E(z) verifies the equation 

1)] E' 



(¥)"«)• 

(13) 



Conversely, (TI3~I) implies the fulfillment of Eq. (JJJ) . 

At first glance, Eq. ( JTBl seems 'much worse' than the original DCHE 
representation. Nevertheless, it is this equation which we shall attempt to 
solve searching for its solution in the form of Laurent series 

oo 

e= a ^ fc ( 14 ) 

'centered' in the points z = and z^ 1 = (which are the only singular points 
for Eq. ([TBI ) with unknown z- independent coefficients a&. Then, carrying out 
straightforward substitution, one gets a sequence of 3-term recurrence rela- 
tions binding triplets of neighboring series coefficients which can be written 



- 3 E" + z[(e- 2ix) z - li (; 



+ 



n — e 



+ ix \ z 



+ (1 - e) 



+ ix 



x 



n + e 

+yU I IX 



E. 



where A = (2a;) — fi . 
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down either as 



or as 



-ft ( k - 1 - - ixj a fc _! 

(n + e \ 
fc + H ixj a fe+1 , (15) 

where = I k H ix 1 — I — - — ] , (16) 



-/x (k - 1 - + ixj a_( fc _i) 

+ (Z k + A)a_ fc + /i ( fc + H h ixj a_( fc+ i), (17) 

where Z k = ( k + + ix) - f (18) 



and = 0, ±1, ±2 .... The sets of Eqs. (fl5l) and (I17p are exactly equivalent 
and each of them covers the whole set of equations the coefficients a k have 
to obey. However, in the approach utilized in the present work, we shall 
consider them both in conjunction, employing ( fl5l) for coefficients a k with 
fc > -1 and flUD for a k with k < 1. Thus, Eqs. (JT5D and (HTJ) will be 
considered separately but on the common index variation 'half-interval' k > 
(remaining legitimate, in principle, for arbitrary integer k). Obviously, 
these two equation sets cover the complete set of conditions imposed to the 
coefficients a k and are 'almost disjoined' intersecting in their 'boundary' k = 
±0-members alone. 

Let us further consider for k > the following formal infinite products 

oo oo 

Rk = l[M j , R k = Y[Mj (19) 

j=k j=k 



of the 2x2 matrices 



M S = gV*,), (20) 
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It is assumed throughout that the matrices Mj, Mj with larger indices j are 
situated in products to the right with respect to ones labeled with lesser 
index values. 

Notice that in the case x = and integer n, zero may appear in denomina- 
tors of M-factors, making the above definitions meaningless. This apparent 
fault admits a simple resolution (see Eq. (1581 and the discussion following it). 
For a while, we temporary leave out consideration of such specific parameter 
choices. 

It is also worth noting that the above definitions of Rk, Rk may be un- 
derstood as a concise form of representation of the 'descending' recurrence 
relations 

Rk = MfcRk+i, (22) 
R k = MkRk+i, k = 0,1,... (23) 

among the neighboring i?'s. These are the only dependencies which will be 
actually used below in derivations involving Rk,Rk- 

The formulas (fl9l) are 'formal' since neither the issue of the convergence 
of such sequences nor how one has to understand the 'initial values' R^, i?^ 
necessary for the actual determination of the 'finite index value' i?-matrices 
are here addressed. 

Now a straightforward calculation applying Eqs. fl22|) . fT23|) shows that the 
following formulas 

a* = /* fc w,, (1 t A + T ix) , (o, i) • * ■ ( I ) (24) 



r (k + i + - ix) 



r(i + ^ + ix ) (1 



r (fc + i + s=£ + ix) 



yield the formal solutions to Eqs. fll5p and (1171) . respectively. 



4 Validation of the formal solution 

In this section we show that the formal solution of Eq. ffl3|) presented in the 
form of expansion (}Tlj) with coefficients given by Eqs. (|24p . (j25p represents a 
well defined analytic function of z. This means, first of all, that the infinite 
matrix products Rk, Rk it involves converge. Moreover, the convergence takes 
place for any constant parameter values. 
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The key auxiliary result which may be utilized for the proof of this asser- 
tion is as follows. 

Lemma. 

Let us consider the sequences of complex numbers atj, (3j,jj,Sj 
satisfying the following 'ascending' matrix recurrence relation 

7, 6 A J ~ I 7,-1 S^x )\ a, ) ^ 



where 



either <jj = 1 

or crj = Zj _ x / Zj . 



(27) 



Then all they converge as j ' ^ oo. Moreover, lim/3,, = = UmSj 
whereas for a, 7-sequences there exist positive quantities N a , iV 7 
and positive integer jo, all depending at most on n, \i, A, x, e, such 
that for all j > jo 



\<*3 



limj'-iooafl < N a m8bx(\aj>\)j , (28) 



?>30 



\jj - lirrif^jfl < iV 7 max(|7 i /|)j 1 , 

j'>jo 

where the maxima are finite. 

The outline of the Lemma proof can be found in the Appendix. 

Remark: Formally, we need not include in the lemma stipulation the re- 
quirement that either x ^ or n is non-integer (which would a priori evade 
possibility of contributions with zero Z* in denominator) because with fixed 
constant parameters and sufficiently large jo no zero Z* may appear. 

Let us return to Eqs. (fT9|) and consider the four double-indexed sets of 
complex numbers {a, (3, 7, 5})f defined as follows: 



7 io) S 



60) )=K =UM k (29) 

k=jo 



It is straightforward to verify that the sequences obtained by the picking out 
the elements with common value of the upper index j obey the recurrence 
relations (126]) for the upper choice in (1271) . Hence it follows from the Lemma 
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that all they converge. We denote the corresponding limits as oft ' etc. We 
have therefore the consistent definition for the infinite matrix products (IT^|) 




(30) 



Let us further note that, increasing j , the 'starting' sequence elements 
( a \i°\ a( j °+i f° r a e ^c) tend to the corresponding elements of the idempo- 
tent matrix 



the discrepancy decreasing as O(j ' 2 ). On the other hand, in accordance 
with (128]) the elements of R^J differ from their j-limits constituting Rj by 
0{Jq ) -order quantities. This means that Rj tends to as jo goes to 
infinity with the difference going to zero as O^ 1 ). In other words, we have 
the 

Lemma corollary: 



Besides, in view of the convergence, 

2. The modules of the elements of all the matrices Rj are 
bounded from above in total 

- provided 'no-zeroes-in denominators' condition for the parameter choice 
is respected, of course. 

In according with the above properties, decomposing Rj into the product 
of the two factors, Rj = Rf -Rj , and approximating Rj by M^, one obtains 
a simple but important algorithm of computation of the products f fT9|) . The 
approximation 




(31) 




oo 



30 



Rj-YlMk-M, 



(32) 



k=j 



is the better, the larger j > j is selected. In the limit, one gets 




(33) 



10 



This interpretation resolves (in a quite obvious way, though) the afore- 
mentioned problem of specification of the 'initial value' for the sequence 
Rj treated through the 'descending' recurrence relation '[Roo] Rj —* 

Rj-i implied by Eq. (1221) . (As a matter of fact, an arbitrary constant 

matrix whose product with M m is nonzero, including the unit matrix, might 
be used instead of in (1521) . affecting only the overall normalization of the 
result and the accuracy of (1521) -like approximation for given 

Having introduced the consistent representation of the matrices (1T9]) . it is 
straightforward to do the same for the matrices Rj (12U1) . The above specula- 
tion applies to them with minor modifications as well. The only distinction 
is the making use of the second (lower) choice in Eq. (12 7p and the operating 
with complex conjugated quantities (equivalent in our case to the replacing 
x by — x) throughout. We shall mark the elements of Rj obtained in this 
way with tildes over the corresponding a's and 7's. 

Now, with the consistent interpretation of the products Rj,Rj ( 119]) in 
hand, one is able to calculate the coefficients a^, for k — 0, 1, 2 . . . in 
accordance with Eqs. ([21]) ([25]). The triple matrix products reduce to separate 
elements of i?'s (or i?'s) denoted above as 7 (*) (for Eq. flU}), and ^(for 
Eq. (1251) . respectively) which are the functions of the parameters n, fi, A, x, e. 
Therefore, the sequences 

^ 7 r(A; + l + ^-ix) 1 J 



k~(k) 



S„ fc = ^P- ~r^- — 2 ' '"' . , fc = 0,l,2,... (35) 

z k _ x r (k + 1 + v + ™) 

are well defined and solve Eqs. (1T51) . (1T7|) . respectively, for k = 1, 2, . . . 

The important feature of the expressions (134|) . (135|) which is used below 
is their asymptotic behavior for large values of the index k which is easy to 
derive in explicit form. Specifically, in accordance with inequalities (T28l) . the 
set of modules of 7- and 7-factors involved in Eqs. (1511) . (1551) is bounded in 
total from above and each of these sequences converge to a finite limit. These 
imply in particular the validity of the estimates 

\a k \ oc ^, |a_ fc | oc (36) 
asymptotically, in the leading order, for some k- independent (1X2- 



r(i + ^ + ix) 
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5 Matching condition and phase-lock crite- 
rion 



By construction, the a-coefficients defined by Eqs. (1241) and fl25|) obey the lin- 
ear homogeneous equations f|T^]) and (1TT1) . respectively, which are 'the same', 
essentially, differing only in the associated intervals of the variation of the 
index, consisting of the positive integers for the former and negative ones for 
the latter. However, these two sequences cannot be joined, automatically, 
since they are 'differently normalized'. This means, in particular, that their 
'edge elements' indexed with zeroes, generally speaking, differ. We will de- 
note them as a Q and a_ , respectively, distinguishing here, in notations, the 
index '0' from the index '—0'. 

Now, referring to Eqs. fl34l) . (1351) . one notes that in view of the factors 
of T-functions present in denominators and leading to asymptotic behaviors 
the following power series in z 



oo ~ 



E + (z) = l + J2 ^> £-(*) = ! + £ (37) 
k=i a ° k=i a '° 

admit absolutely converging majorants. (We assumed above a± ^ 0. Oth- 
erwise, i.e. if a±o = 0, a±\ may not vanish and the series with the coefficients 
a±fc/a±i can be utilized instead.) Indeed, the Maclaurin series for the expo- 
nent can play this role. Therefore, the series E + (z) and -E_(;? _1 ) define some 
entire functions of z. As a consequence, the expression 

= ig L a^M^g(^M + E _ M _ d 

(n + e — 2ixJ (n + 2 — e + 2ixJ 

(38) 

represents a single- valued function analytic everywhere on the Riemann sphere 
except the poles z — 0, z~ x = 0. They are the essential singular points for E. 

It has to be noted that the additional ^-independent fractional factor 
in (1381) may be regarded as a specific common 'normalization' of the (I37p - 
type series which may be, in principle, arbitrary. However, its given form is, 
essentially, unique as being fixed (up to an insignificant nonzero numerical 
factor) in view of the following reasons. 

The two sine-factors in the numerator regarded as holomorphic functions 
of n + e ± 2ix are introduced for the canceling out zeroes in denominators 
arising due to the poles in the factors Z~ l and Zj 1 involved in the products 
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and regarded as the functions of the same parameters. The set of these 
(vicious, essentially) singularities constitute a homogeneous grid which is just 
covered by the grid of roots of the sine-factors in fl38l) — with the two ex- 
ceptions. These two 'superfluous' sine-factor roots are, in turn, 'neutralized' 
by the two linear factors in the denominator in fl38l) which are therefore also 
uniquely determined. As the result, in vicinity of any zero in denominators 
in coefficients of the power series defining E(z) considered as the function of 
n + e ± 2ix (a root of some Z* or Z\ ), the resulting expression takes the 
form of the ratio sinx/x (x ~ 0) and is not now associated with any irregular 
behavior. Thus, as a matter of fact, the fractional factor involved in (I3"8j) is 
distinguished (up to a numerical factor) by the claims (i) to cancel out the 
poles in the original expressions of the a-coefficients (13"4"j) . (|33|) considered as 
the functions of n + e±2ix and (ii) to introduce neither more roots nor more 
poles as a result of such a 'renormalization'. 

Now, when plugging the function fl38l) in Eq. fflBl in order to verify its 
fulfillment, we may provisionally drop out, sparing the space, z- independent 
fractional factor (restoring it afterwards). 

It is important to emphasize that the expressions (1341) . (1331) . by con- 
struction, verify all the 3-term recurrence relations (fl5l) . f|T7|) which bind the 
a-coefficients with indices of a common sign, either non-negative or non- 
positive. The only equation which does not fall into the above categories, 
and, accordingly, has not been automatically fulfilled, is the 'central' one 
binding the coefficients e. the equation 

/ 77, — 6 \ f Tl ~\~ € \ 

H 1 + —r- + ix a_i + (Z + A)a + n 1 + — he \a x = 0. (39) 



With normalization adopted in fl37|) . one has a = 1, ai = ai/a ,a_i = 
a_x/a_o- Further, in accordance with (J31|) . ([35]) 

Besides, one has = a^ 1 ^, = 6ft*'Z_i/ Zq. Combining these dependen- 
cies, the following representation of Eq. f[3"9"j) arises 

= ^ + (Zo + X)+^. (40) 
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Accordingly, it is convenient to introduce the following function of the pa- 
rameters x and n, A, //, e 

4r sin (|(n + e - 2ix)) sin (|(n + e + 2ix)) 

(n + e - 2ix) (n + 2 - e + 2ix) 
(/i 2 7 (1) « (1) + (Z + A)5«a« + ^W a W) (41) 

where Z = ((e - l)/2 - ix) -((n + l)/2f (see (116]) ) and a's, 7's are defined 
as the elements of the convergent matrix products as follows: 

(#S)-nM$S)-fl* (42) 

(see Eqs. ( |T9l) -(l211). The fractional multiplier in the first line of Eq. ( ]4"Ti) 
coincides with the one entering Eq. (|38|) and plays the identical role: it 
eliminates the vicious singularities arising for specific values of the parameters 
n, x. We shall name H = S(x, n, /i, A, x, e) = H(x;...) the discriminant 
function for brevity. The following statement holds true. 

Proposition 

Restricting x to real values, the equality H(x; . . . ) = is the 
necessary and sufficient condition for the single valued analytic 
function (|38|) to verify Eq. (fT3l) everywhere on the Riemann sphere 
except its poles 2: = 0, z^ 1 = 0. 

Indeed, the vanishing of H implies the fulfillment of ( |39i) (where a's are 
expressed through a's), the last equation binding coefficients of the expan- 
sion (|T4|) which has not been fulfilled as the result of the very coefficients 
definition. Now all the 3-term recurrence relations for a-coefficients, which 
Eq. (fT3|) is equivalent to, are satisfied and the analytic function (1381) verifies 
Eq. fll3p everywhere except of its own singular points z = and z~ l = 0. 

The equation 

S(x;...) = (43) 

referred to in the above proposition can be named the matching condition 
since it enforces the sequences of the coefficients a^, a_fc, separately obeying 
the corresponding 'halves' of the equation chain f|T5l) (equivalently, ffTTj) ) to 
be 'matched' in their 'edge' elements a ±0 = l,a±i- 
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It is worth emphasizing that, up to this point, the parameter x (absent 
in Eq. (J7J) has not been restricted in any way (it was only assumed to be 
real). Now and in what follows we regard the condition (j4"51) as x definition 
eliminating this odd 'degree of freedom'. Now it is a well defined function 
of the other parameters. It seems interesting enough that the addition of 
unspecified constant x to the transformation f[T2"l) and its subsequent 'fine 
tuning' by means of the claim of fulfillment of Eq. fl4"3"j) is necessary for the 
representation of solution of DCHE (J7J) in terms of convergent Laurent series. 
More precisely, it is clear that Eqs. (fl~6l) can be solved for any (including 
trivial zero) choice of x, choosing loosely cij , a JO +i for arbitrary jo an d then 
calculating, term by term, all the coefficients aj, advancing in parallel in both 
directions of j-index variation 'from j towards ±00'. Then (I14p immediately 
yields a (x-dependent!) formal solution to Eq. ffl~3l) and hence, through 
transformation ([12]) . to Eq. ([7]). However, it can be only formal and will 
necessarily diverge for any z unless the matching condition (H3j) is fulfilled 
- just in view of the uniqueness of solution with the analytic properties 
presupposed. On the other hand, considering separately the 'halves' of the 
set of Eqs. ( 1161) and resolving them 'in index variation directions' opposite 
to the ones assumed above (in a sense, 'from ±00 towards ±0'), we obtain 
the always converging series (|37|) . However, as we have seen, we again have 
no solution (in this case even formal) unless the matching condition fixing x 
is fulfilled. Obviously, the uniqueness property implies that the introduction 
of the 'branching' power function factor, as in (1121) . is the only way to obtain 
a solution to Eq. ([7]) admitting representation in terms of convergent power 
series. 

Now, tracking back the relationship connecting Eq. (fl~3]) with the pri- 
mary Eq. ([T]) and invoking the general theory of the latter applicable in the 
case of arbitrary continuous periodic q(t) one can infer the following 
statements which however, in the present context, are only of the status of 
conjectures in view of the lack of their proof 'from the first principles', i.e. 
on the base of the properties of the discriminant function S following from 
its definition. 

Conjecture A. 

1. There exists an open non-empty subset of the space of the 
problem parameters n, fi, uj, e where Eq. (1431) admits a real valued 
positive solution. 

2. If real x solves Eq. fl4"3"l) . — x is the solution as well. No more 
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real roots exist. 

3. Real roots of Eq. (T4"3"i) obey the condition |x| < (2u)~ 1 . 

Remark: The last statement is nothing else but the form of the limitation on 
the rate of growth or decreasing of the functions x, y = x, y\ z=e ^t of the real 
variable t implied by the inequalities ([6]) and the equation 



x — iy = —iz 




following from definitions. Notice that the latter clarifies the role of the 
parity parameter e which determines the multiplicity of the inverse to the 
map S 1 3 z — > [x + iy)/\x + iy\ G S 1 induced by the solution (I3"8|) . If 
e = 0, the revolution along the circle \z\ = 1 leads to the reversing of the 
direction of the vector with components (x, y) whereas in the case e = 1 its 
direction is preserved. The inverse map is double-valued in the former case 
and one-to-one in the latter one. 

More properties of the discriminant function can be inferred from the nu- 
merical experiments although, as opposed to the assertions of the Conjecture 
A, they have, to date, no analytic arguments in their favor yet, even indirect. 
Nevertheless, the first item below is important enough from viewpoint of ap- 
plications (seeming also plausible enough) to be explicitly formulated here. 

Conjecture B. 

1. Phase-lock criterion. 

Equation (|43[) admits a real non-zero solution if and only if 

3(0;...) >0. (44) 

This means in particular that 3(0; . . . ) is real; moreover, the numerical study 
makes evidence that 

2. 3(x; . . . ) is real for real x 

(assuming the other parameters to be also real, of course). 

6 Floquet solutions of DCHE and involutive 
solution maps 

Let us assume now that there exists a real positive solution x of Eq. fT43|) . 
With this x, the function E(z) defined by Eq. (1551) verifies Eq. ffT3"j) . Let us 
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consider the function E#(z) defined through E(z) as follows: 



E. 



„2ix~e 



Z 



+ 



n + e 



ix I - // ) E ( - 



(45) 



Then a straightforward computation shows that it verifies Eq. ( |T3l) . provided 
E(z) does. 

It is worthwhile to note that the repetition of the transformation (1451) 
yields no more solutions to Eq. (TT3T) . As a matter of fact, one has # o # = 
{2ui)~ 2 Id. Thus (2c<j) _1 # is the involutive map on the space of its solutions. 

Next, the functions E(z) and E#(z) are linearly independent for nonzero 
x. Indeed, utilizing ffT3l) . one finds 



v 2i»r-e-l 



n + e 



ix + fiz ) E' [ - J + 



+ I I — i ix 



> 2 + 1) + 



n 



— IX 



+ A \z 



El 1 - 

z 



.(46) 



This reduction allows one to calculate the determinant of the linear trans- 
formation binding the pairs of functions E#,E#'(z) and E,E'(l/z) which 
proves equal to (2uj)~ 2 : z 2 (- e + 2lJ< ) anc } j s therefore nonzero. Hence E#(z) is 
not identically zero (and may vanish at isolated points, at most, as well as 
E(z)). Finally, E(z) is periodic on the unit circle centered at zero whereas 
E#(z), for real x ^ 0, is not. Hence they are linear independent. The 
functions E(z) and E#(z) constitute therefore the fundamental system for 
Eq. (fT3|) and any its solution can be expanded in this basis with constant 
expansion coefficients. The analytic properties of these solutions identify 
them as the unique pair of the Floquet solutions of the reduced DCHE under 
consideration. See [T4] , section 2.4. 

The function E(z) obeys the important functional equation which can be 
derived as follows. A straightforward calculation shows that the RHS expres- 
sion of Eq. (|45l) with the 'branched' factor z 2l>c removed (i.e. z~ 2l *E^z)) 
satisfies the ODE which coincides with (TT5|) up to the opposite sign of the 
parameter x. This means that, for real x, n,fi,u, the analytic function 



E(z) 



E'[- z ) + 



n + e 



+ ix ) z 



(47) 



where E(z) = E(z), is the solution of Eq. ( |T3|) itself. 
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As opposed to E#, this 'yet another' solution E(z) has the same analytic 
properties as E(z) and hence must coincide with it up to some numerical 
factor Cc, i.e. 



C c E{z) 



(48) 



(Cc may not vanish since otherwise E# would also be zero.) This is the 
generalization of the similar property of the so called 'Heun polynomials' 
established in Ref. [T2] . 

The complex valued constant Cc actually reduces to a single real con- 
stant. To show that, let us notice at first that if E(z) verifies Eq. (|T3l) it 
follows from the latter and 



C c E\z) = z 



-e-1 



" ( ix + fiz) E\l/z) 



2 



+ [ fi e -±^(z 2 + l)-ifix{z 2 -l) 
+ ( - - 2 + A) z) E(1/z) 



(49) 



Evaluating now Eqs. (1481) . (1491) together with their complex conjugated ver- 
sions with z = z^ 1 = 1, one obtains four linear homogeneous equations 
binding the quantities E(l), E'(l), E(l), E'(l) which may not vanish simul- 
taneously. The corresponding consistency condition reads \Cc\ 2 = (2co>)~ 2 
implying 

C c = (2co)- 1 e iC % (50) 

where C c is some real constant (actually, the function of the parameters 
n,(/,,w,e). It encodes all the monodromy data for Eq. (fT3|) . essentially. 

It is straightforward to show that the transformation Eq. (l4*8l) is also 
involutive. It manifests the specific symmetry in the behaviors of the function 
E(z) in vicinities of the essentially singular points z = and z^ 1 = 0. 
Remarkably, this symmetry implies itself the fulfillment of Eq. (1131) . Indeed, 
differentiating (Jl8"j) and taking into account (150]) . one arrives at Eq. (|T3l) . In 
a sense, Eq. ( |4"8]) together with stipulation for the analyticity of E(z) can 
be considered as the equivalent to Eq. (|T3|) . Additionally, Eq. (jlSj) implies 
anti-linear (involving complex conjugation) dependencies among the 'distant' 
Laurent series coefficients and 0^,(1^+1 (whereas Eq. ([TBI) (or (|17[) ) binds 
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'nearby' ak^k^i)- m particular, it suffices to find all for k > and then 
a_fc can be computed from the latter by means of a simple transformation. 



7 Essentially periodic and general solutions 
of overdamped Josephson junction equa- 
tions 

The connection between the functions E(z), E#(z), E(z) pointed out above 
is important for the lifting the results concerning solutions of Eq. f[T5j) to the 
level of original OJJE. This procedure applies Eqs. (JH), flSJ), (fTUj) . ffT2"j) and 
leads to the following conclusions. 

At first, the representation of the two special (and the most important) 
solutions to OJJE for which the exponents exp(i</?) are periodic (for brevity, 
we shall call such phase functions essentially periodic) follows. It reads 

( E '(z) n + e \ 
exp(-i^) = 2luJ [ z -j^j + — 2 (51) 

exp{vp) = -2iu Iz- 1 E K ' + -ix -nz J , (52) 

where z = exp(ia;t). 

For x > the first of these formulas determines the asymptotic limit (the 
attractor) of a generic solution whereas the second solution is unstable (the 
repeller). It is important to emphasize that the functions tp(t) defined by 
Eqs. ( T5Tj) and (!52|) are real and Eq. (!50j) is the crucial property utilized in 
the calculation establishing this fact. 

At second, it is straightforward to obtain the 'nonlinear superposition' of 
solutions (f5"Tj) . (1521) operating with their DCHE-related counterparts. The 



19 



result is represented by the formula 



exp(i^) 



- < cos'0 • E(z) + smip ■ z 



E\z- 1 ) + ( ( - ix ) z - fi ) E(z~ v 



X 



UJ cos ib ■ 



zE'(z) + ( — ix-nz) E(z) 



H sin ib ■ z 



-1 



(53) 



where ib is an arbitrary real constant. More exactly, the set of all functions 
ip described by Eq. (1531) is parameterized by a point on the unit circle. As 
opposed to fl5T|) . fl52l) . the function ff^9l) is defined on the universal covering 
of the Riemann sphere with punctured poles, Q. Continuous (and then real 
analytic) function ip(t) determined by this equation on Qi G £1, where z is 
understood as e luJt , is just the general solution of OJJE in the case of phase- 
lock. 

In particular, Eqs. (!5TT) . (l52l arise as particular cases of (|53|) for ib = 
and ib = tt/2, respectively. As a consequence, asymptotic properties of the 
general solution mentioned above immediately follow. Indeed, as t increases, 
the exponent (1531 is converging to (15T|) and is moving off (1521) (unless it 
coincides with the latter). The two solutions described by fl5TT) . fl52l) are the 
only ones which are not affected by the translations t — > t + 2-koj~ 1 (in the 
sense the exponents (I51l) . fl52l are kept unchanged) and preserve their form 
in asymptotics. 

At third, considering ip defined by fl5"T|) as analytic function of z and 
taking in account Eq. (jTBl . one obtains 



dip 
dz 



-iz 



- I z + z 

E'(z) 
z- 



n 



E z) 



\K I — \l 



n + e 



- ix J — fiz + 



4u 2 



x 



\IZ 



(54) 



On the unit circle, this p> verifies OJJE. It is therefore smooth (even real ana- 
lytic). Then (I5"4"|) is continuous on the unit circle. Finally, since exp i<£>| 2=e xpiurf 
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is periodic, the following proposition holds true: 

Proposition 2. 

The quantity 

k = (2ir)- 1 j^-dz, (55) 

where dip/dz denotes the RHS expression from Eq. (j5"lj) and the 
integration is carried out over the circle \z\ = 1, is well defined 
and equals to an integer. 

This integer is the degree of the map S 1 S 1 induced by the function fl5Tj) . 
In physical applications, it is called the phase-lock order and is considered as 
an integer- valued function of the parameters. Phase-lock order is involved in 
the formula representing the property of being 'essentially periodic' for the 
phase function defined by Eq. (|5T|) (and asymptotically for a generic phase 
function) which reads 

Vt (p(t + 2-KOJ- 1 ) = <p(t) + 2nk 

In a phase-lock state of JJ, the uniformly distributed discrete levels of aver- 
aged voltage equal to k ■ (hw/2e) for some k = 0, ±1, ±2 . . . are observed. 

Conjecture C. 

Any integer map degree (155]) is realized on some non-empty open 
subset of the space of the problem parameters n, /i, u>, e. 

This assertion is closely cognate to the item 1 of the above Conjecture A. 

8 Summary 

It the present work, the general solution of the overdamped Josephson equa- 
tion (JTJ is derived for the (co)sinusoidal RHS function (j2J) in the case of one 
of three possible asymptotic behaviors known as the phase-lock mode. The 
solution is represented in explicit form in terms of the Floquet solution of 
the particular instance (corresponding to the vanishing of one of the four free 
constant parameters) of the double confluent Heun equation (DCHE). The 
Floquet solution of DCHE is represented in terms of the Laurent series whose 
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coefficients are determined by the convergent infinite products of 2 x 2 ma- 
trices with a single zero element tending to the idempotent matrix (13 ip . The 
derivation presupposes the existence of a real solution of the transcendental 
equation (|43p which is equivalent to the claim of realization of the phase-lock 
mode for the given parameter values. The plausible criterion of its existence 
(i.e. the phase-lock criterion) is conjectured. 

It is worth summarizing here the main steps of solution of OJJE. They 
can be condensed as follows. 

• The investigation of the basic properties of Eq. ([I]) for arbitrary peri- 
odic (sufficiently regular) q(t) allows one to establish the division of the 
space of the problem parameters into the two open areas of which one 
corresponds to the phase-lock property of the OJJE solutions whereas 
another corresponds to their pseudo-chaotic behavior revealing no sta- 
ble periodicity. For the (lower-dimensional) complement to these areas 
the intermediate behavior is observed. The corresponding results are 
discussed in sufficient details in Ref. 

• The next important point is the intimate connection (first mentioned 
by V. Buchstaber, see, e.g., Ref. [T3]) between (jTJ and a simple linear 
system of the two first order ODEs flSJ). For (co) sinusoidal RHS function 
([2]), the latter takes the form 

• At the next step, the transformation was found which converts the 
linear system to a particular instance of the double confluent Heun 
equation (TTTj) . 

Generally speaking, it could be solved by means of the expansion in Laurent 
series [2], [12] centered at the singular points but preliminarily the additional 
simple but important transformation has to be carried out: 

• the 'branched' power factor involving unspecified constant (x-dependent 
contribution in Eq. f|T2l ) is introduced. 

• The addition of the discrete 'parity' parameter e, assuming either the 
value or the value 1, which is involved in the power factor proves 
necessary for the subsequent exhaustive 'indexing' of the solution space. 

• After that, the standard technique of the power expansion leads to the 
'endless' sequence of the 3-term constraints (fT5l) (or, which is the same, 
( HTTP ) imposed on the unknown series coefficients. 
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• The next step is the devision of the set of power series coefficients 
into two subsets. The non-negative-index- value coefficients and non- 
positive-index-value ones are treated separately, solving the separate 
subsets of the equations (fl5j) and ffTTj) . respectively, for k > 1. The 
application of the continued fraction technique leads, after some trans- 
formations, to the 'explicit' formulas for the series coefficients involving 
infinite products of 2 x 2 matrices converging for large index values to 
the idempotent matrix (I3TT) . This convergence is sufficiently fast to 
imply the convergence of the matrix products and, accordingly, the 
finiteness of the series coefficients. Moreover, the associated estimates 
make evident the existence of the absolutely converging majorants for 
the resulting Laurent series. Therefore, they actually determine the 
Floquet solution of DCHE. The latter proves representable as the sum 
of the two entire functions of the arguments z and z~ x , respectively. 

• The procedure producing Laurent series coefficients noted above proves 
suffering however from the improper introduction of a kind of vicious 
singularities arising as zeroes in denominator which appear for some 
special parameter values. They are eliminated my means of multipli- 
cation of the 'raw' coefficient expressions by some z-independent (but 
parameter dependent) factors given in explicit form. 

• Now the 'solution candidate' for Eq. ffl3l) can be represented as the 
analytic function (|3"8|) which is well-defined for any parameter values. 
However, at the price of automatic convergence of the series it has been 
built upon from, it does not always verify Eq. (fi"3]) . The equation is 
fulfilled if and only if the fulfillment of Eq. fH3|) , which is the transcen- 
dental equation for the still unspecified parameter x, is stipulated. 

At that stage, having solved Eq. (143]) . a single solution (the Floquet solution) 
of DCHE can be regarded as having been explicitly constructed. 

• The invariance of the space of solutions of DCHE under consideration 
with respect to transformation represented by Eq. (f45|) allows one to 
immediately obtain the fundamental system of its solutions in terms of 
the single Floquet solution noted above. 

• The automorphism represented by Eq. (HHI) expresses the important 
intrinsic property of the Floquet solution of DCHE. It is used for the 
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derivation of the explicit representation of the exponent exp(i</?) speci- 
fying the real valued phase function ip to be obtained as the restriction 
of the analytic function from the universal covering of the Riemann 
sphere with punctured poles (Eq. fl53|) ) to the universal covering of the 
unit circle. It yields the general solution of Eqs. ([I]),© in the case of 
phase-lock. 

• Employing analytic properties of exp(i(p)(z), the formula fl55|) involving 
Floquet solution of DCHE follows which gives the degree of the map 
S 1 =^ S 1 it induces (the winding number) also known in application 
fields as the phase-lock order. 

It is worth noting in conclusion that all the constructions derived above 
admit a straightforward algorithmic implementation which have been used 
for the numeric verification of the relevant relationships. 
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A Appendix: outline of the Lemma proof 

Eq. (126]) implies (3j = i^aj^Zj ,Sj = fi 2 'jj_ 1 Z^ 1 and hence the asserted 
properties of the sequences /3j,5j follow from the existence of finite limits 
for the sequences atj,jj. As to ctj and jj, they have to obey the identical 
decoupled 3-term recurrence relations which, for a's, read 

ay = (1 + \Zjl {e _ l)j2 )a^i + ^ 2 Zr 1 . +(e _ 1)/2 a i _ 2 (56) 
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where e = 1 for the upper choice of a in (j27|) and e = for the lower a choice 
therein. It suffices to consider the a-sequence case. 

Evidently, for every integer jo > and I > any solution of equations 
( 156]) can be represented in terms of the decomposition 

a h+i = C 1 + Pj ,i)ocjo-i + qjo,iocj -2, (57) 

for some coefficients Pjj, qjj independent on the 'starting' terms azj -i, o(j -2- 
Applying fl56l) . it is straightforward to show 

Qh,i+i = ( 1 +Ph+t,i)f j2z io-^ ( 58 ) 
whereas for p*^ one gets the following recurrence relation: 

Pjo,i+i = Pjo+i,i + M 1 + Pj +i,i) z J X + A 1 + Pio+2^-i)^ r +i- e - (59) 

Eq. fl59|) is equivalent to fl56|) . essentially, but it possesses the advantage of 
being endowed with the standard 'initial conditions' 

Pio-i = 0,p iOi _ 2 = -l (60) 

which follow from definitions. Besides, one gets 

9io.-i = 0,g iOi _ 2 = l. (61) 

It proves convenient to carry out one more rearrangement of unknowns 
introducing the differences 

APio.i = Pj ,i+i ~ Pj ,i- (62) 
which obey the own 'initial conditions' 

Ap i0l - 2 = 1, (63) 
Ap,o,-i = AZ" 1 , (64) 

and similar recurrence relations 

A Pj ,i+i = A Pjo+i,i + ^Apjo+yZT 1 + /i^^+a^!^^.-. (65) 
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Now, summing up the subset of the latter with the common sum of indices 
at the left and taking into account (1641) . all but two 'free' Ap-terms cancel 
out and one obtains the following equation 

l+l 



+ A ^2 APj +l+m,Z-mZ JO + m 
m=0 
l+l 

(66) 



m=0 



In the sums, the second index of Ap* i# is everywhere less than the same index 
at the left that allows to apply the method of mathematical induction. For 
the 'starting' values -1,0 of the second index one has 

Z jo Ap jo ,-i = a, 

z k+i^Pkfi = A(l + AZ^ 1 ) + n 2 Zj 0+1 Zj Q+2 _~. 

Therefore for I = —1,0 there exist the finite limits Hindoo \Z, +l+1 Apj 0) i\. 
As a consequence, for these Vs one has 

(67) 

for appropriate constant iV which is convenient to choose > 1. Let us consider 
this fact as the starting point of mathematical induction and assume that 
for some integer l > and any integer I from the interval [— 1 , Iq] , (1671) 
holds true. We may apply it for the estimating from above of the quantity 
\Apj (h i 0+ i\. This can be realized making use of the 'decomposition' fl66l) and 
the following elementary inequalities 

L+jo+l 1 _)_ U_|_ ]|-1 

2 |Zm|_1 < Z _lr,-Lll/9 x/ f -1)/9' (68) 



m=J 
L+jo+2-e 



m=jo+l-e 



Jo -|n + l|/2 + (e-l)/2' 

1 + \n + 1| _1 
j -|rc + l|/2 + (e-l)/2' 



T \Z m \- 1 < + w , (69) 
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where L > (and -1) . These imply the inequalities 

\Ap j0th+1 \ < \X\\Zj, 



h+1 



J jo+m+l-e\ 





(|A| + |/i 2 |)(l + 


n + 1| 


- 1 ) 


l^j'o+io+ll 


(jo - 


n 


+ l|/2 + (e-l)/2) 



+iv|^ Wo+1 |- 1 iai ^ i^r 1 + i^ 2 i 

\ m=0 m=0 

,Z,.,,.,ol (IAI 4- l// 2 l 

Since we assumed N > 1, there exists the lower index value bound such that 
for any jo exceeding it the factor in brackets is less than N and then the above 
inequalities imply \Apj 0t i 0+ i\ < A r |Zj 0+ / 0+2 | _1 . (|67j) is therefore established 
for sufficiently large jo an d arbitrary I > 0. Increasing iV if necessary, f l67l) 
proves valid for arbitrary jo- 
in view of this property, one sees that the sum Y^lo APio,J nas the ma j°- 
rant ^ z |Zj 0+ / +1 | _1 and thus converges itself. The sequence of its partial sums 

EL=o APjo,m = Pj ,i+i -Pj ,-i = Pio.J+1 also converges as / -> oo. Moreover, 
in view of ([SO]) . ([52]) . ([55]) . (I6TI) one has the Z- uniform bound 

r r v^, , 1 iV(l + \n + lh 1 ) , . 

l ft „ J+1 | < w £ IW < ^TTT(^ip - (70) 

m=0 J v " 

It follows from Eqs. (I57j),(l58l) 

a i = a h+i = ( 1 + PooMh-i + (1 + Pj +i,l-i)^ Z i -e a h-2 (71) 

and the convergence of a-sequence follows from the convergence of p* t i as 
I — > oo. Then one has 

lim otj- - a io _i = limp i0jZ a io _i + (1 + limp io+ i,j)/?Z^a jo _ 2 - (72) 

The factors in front of the first and second terms to the right scales as j^ 1 
and jo~ 2 , respectively. We may therefore write down the following inequality 

| limaj - a jo -i\ = jVmax(|aj _i|, l^o-aDj'o" 1 * 

where N may depend on the parameters n, A, /x, x, e but not on the specific 
specimen of a-sequence. This obviously implies the inequality (1281) . 

It has also to be noted in conclusion that the case n = — 1 formally falling 
off the above speculation does not actually correspond to an exceptional 
situation. Although inequalities (|68p formally fail, similar ones differing from 
(I68p in the values of 'constant' (jo-independent) terms alone can be derived. 
The further reasoning holds true and leads to the same conclusions. 
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